Pion interferometry for hydrodynamical expanding source with a finite baryon density 
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We calculate the two-pion correlation function for an expanding hadron source with a finite 
baryon density. The space-time evolution of the source is described by relativistic hydrodynamics 
and the HBT radius is extracted after effects of collective expansion and multiple scattering on 
the HBT interferometry have been taken into account, using quantum probability amplitudes in 
a path-integral formalism. We find that this radius is substantially smaller than the HBT radius 
extracted from the freeze-out configuration. 



PACS numbers: 25.75.-q, 25.75.Gz 



The Bose-Einstein correlation of identical bosons pro- 
duced in high-energy heavy-ion collisions, also known as 
the HBT effect, is an important tool for the study of the 
space-time structure of the emitting source As the 
source expands, cools, and freezes out, it is important to 
know what source distribution the HBT interferometry 
measures. Is it the freeze-out source distribution, the ini- 
tial source distribution, or the initial source distribution 
modified by absorption and expansion? The conventional 
viewpoint is that the HBT interferometry measures the 
freeze-out configuration because rescattering of source 
particles are assumed to lead to a chaotic source. How- 
ever, the extracted experimental HBT radii are insensi- 
tive to the collision energy, and i? ou t /-Rside ~ 0.9 — 1.1 
at RHIC 0, • These general features cannot be under- 
stood within the conventional viewpoint 3 0, 01 ■ 

The validity of the conventional assumption on HBT 
is recently questioned as it was pointed out that because 
HBT is purely a quantum-mechanical phenomenon, the 
effects of rescattering and collective dynamics must be in- 
vestigated within a quantum-mechanical context Q- As 
fully quantum treatment of heavy- ion collisions is beyond 
the scope of our present investigation, we shall describe 
the gross dynamics of the hadron system with a finite 
baryon density by classical hydrodynamics. In this let- 
ter, we follow Wong's work Q and take into account the 
effects of collective expansion and multiple scattering by 
constructing the quantum probability amplitude using 
Glauber's multiple scattering model and the trajectories 
of test particles in the path-integral framework. 

The hydrodynamics and the composition of the fluid 
depend on the collision energy. We study in this work 
the dynamics of heavy-ion collisions in the energy range 
corresponding to about 10 GeV per nucleon beam energy 
in fixed target collisions in which nucleons and deltas play 
an important role and pions are produced dominantly 
from deltas and absorbed by nucleons. 

Relativistic hydrodynamics has been extensively ap- 
plied to high-energy heavy- ion collisions Q . The energy 
momentum tensor of a thermalized fluid cell in the ccntcr- 
of-mass frame is 

T» u (x) = [e[x) + p(x)] u»(x)u v {x) - p(x)g fiU , (1) 



where x is the space-time coordinate, e, p, and it M = 
7(1, v) are respectively the energy density, pressure, and 
4-velocity of the cell, and is the metric tensor. The 
local conservation of energy and momentum can be ex- 
pressed by 

^T^(x)=0, (1/ = 0,1, 2, 3). (2) 
The conservation of baryon number gives 



d^f(x) = 0, 



(3) 



where j M = nb(x)u^ is the four-current-density of baryon 
(rib is baryon density). In the thermalized cell, the den- 
sity number ni of the particle species i, the energy den- 
sity £j, the pressure pi, and the local velocity of sound c s 
can be obtained as a function of local temperature T(x) 
and local chemical potential /i^(x) by assuming an ideal 
hadron gas. For simplicity, we consider a hadronic gas 
consisting of nucleons, A(1232), and pions with spherical 
geometry as in Ref. @|. For spherical geometry Eqs. © 
and J3J become 



dtE + d r [(E + p)v] 
d t M + d r {Mv + p) 
d t N + d r (Nv) 



-F, 
-G, 
-U, 



(4) 
(5) 
(6) 



where E = T 00 , M = T or , N 



F=— (E + p), G=—M, U=—N. (7) 
r r r 

In order to bring out the salient features of the effects 
under consideration, we shall carry out model calcula- 
tions for some idealized space-time configurations. In 
these model calculations, we assume that the initial ve- 
locity v(0,r) = 0, and the initial energy density is a 
Gaussian distribution e(0,r) = eo e~ r / 2R ° ^(J- I n a 
local equilibrium cell, we have fj, n = 0, lin = MA = 
iib- Using the HLLE scheme || and with the relations 
of e(T(x),/j, B (x)), p(T(x),Hb(x)), n b (T(x), hb(x)), and 
c s (r(x),/is(x)), we can get the solution of the hydrody- 
namical equations for F = G = U = 0. Then, using the 
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FIG. 1: Temperature (a) and velocity (b) profiles for the 
expanding source at t n = 3n fm/c. 



Sod's operator splitting method we can obtain the 
solution for Eqs. 0,0, and JHJ. 

Figure 1 shows the temperature and velocity for a 
spherical expanding source with eo = 0.5 GeV/fm 3 and 
Rq = 4.0 fm. The corresponding initial temperature and 
baryon density at r = are T = 137 MeV and ribo — 0.36 
fm" 3 , respectively. They are obtained by using the state 
of equation of mixed ideal gas and eo. The tempera- 
ture To = 137 MeV is lower than the estimated critical 
temperature, T c = 160 MeV, of the transition between 
quark-gluon plasma and hadronic gas [9|. 

The two-particle Bose-Einstein correlation function is 
defined as the ratio of the two-particle momentum dis- 
tribution P(ki,k2) to the product of the single-particle 



momentum distribution r{ki)r{k 2 ). ror an expanding 
source, P(k) and P{k\, fc 2 ) can be expressed as Q 

P(jfe) = j d 4 xe- 2Im ^^ x) p{x)A 2 ( K (x) 1 x) , (8) 

P(k ll k 2 )=Jd 4 X 1 d 4 X 2 e -2W.(*i) e -2Z»i *.(*») 

x p(xi)p(x 2 )|$(«i«2 = xix 2 -> x d ix d2 )\ 2 , (9) 

where A(k(x),x) is the magnitude of the amplitude for 
producing a pion with momentum k at x, which is pro- 
portional to the Bose-Einstein distribution characterized 
by the local temperature in the local frame of the cell 



at x. p(x) is the pion-source density, which includes a 
primary source and the secondary source from A decay, 



p(x) 



prim 



(x) 



.second 



Or) 



n 7r (x)S(t) +TriA(x) , 



(10) 



where T = 120 MeV is the width of A. e -2Xm^.(x) is thc 
absorption factor due to multiple scattering 0, EI > 



<J a hs{^/s 1TN )n N (x)dl , (11) 



e _aim *' ( *' =exp 



where cr a bs(v 's-kn) is the absorption cross section of 
7r + N — > A at the center-of-mass energy ^/s^n and dl is 
the spatial line element along the path of particle propa- 
gation. In Eq. @, $(ki«2 : x\x 2 — > XdiXd 2 ) is the wave 
function for two-particles produced at X\ and x 2 with 
momenta K\(xi) and k 2 {x 2 ), and detected at Xdi or Xd 2 
with momenta fci and fe 2 , respectively, 

<J>(kik 2 : x\X 2 -> x d \x d2 ) 
1 



-j=-{A(ki(xi), xi)A(n 2 (x 2 )x 2 ) 



x exp 



i I Ki(x') ■ dx' — ik\ ■ (xdi 
i I k 2 (x') ■ dx' — ik 2 ■ (xd 2 - xj 2 ) 



x exp 

+ A(k 1 (x 2 ),x 2 )A(k 2 (xi),xi) 
x exp 



x exp 



Ki(x') ■ dx' — iki ■ (xdi — x'f 2 ) 
k 2 (x') ■ dx' - ik 2 ■ (xd 2 - x'fi) 



(12) 



where xt\ and xj 2 are the freeze-out points correspond- 
ing to the particles 1 and 2 produced at x\ and x 2 and 
detected at Xdi and Xd 2 , respectively. x'f X and x'j 2 are 
the freeze-out points corresponding to the particles 1 
and 2 produced at x\ and x 2 and detected at Xd 2 and 
Xdi, respectively. As we consider mainly two particles 
whose momenta are nearly parallel, we approximately 
have x'f! — xji and x'j 2 = ip @ 0|. From Eqs. 
©, ©, and lO the correlation function C(fci,fc 2 ) = 
P(ki,k 2 )/P(ki)P(k 2 ) can be written as 



C(fci,fc 2 ) 



d 4 X g i ( fe i _fe 2)-a:+i0e(2:,fcife2) 
2 



xe 



— 2 Im tfisiz) 



p e s{x; ki_k 2 ) 



(13) 



where p e ff is the effective density 



Pea (a:; feife) = 7==^=== , (14) 



v /P(fc 1 )P(fe 2 ) 
/init(K(a;),a;) = ,o(x)A 2 (k(x), .x) 



(15) 
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FIG. 2: Two-pion correlation functions C(q) for the cases 
(a) freeze-out, (b) present formulation including absorption 
of pions by nucleons, and (c) without absorption. 



is the phase space distribution of particle sources, and 



(f> c (x, kik 2 ) 



{[ Kl {x')- K 2 (x')}-[k 1 -k 2 }}-dx'. 

(16) 



It can be seen that the two-pion correlation function is 
related to the phase space distribution of the pion pro- 
duction source, modified by an absorption factor arising 
from multiple scattering and a phase factor cj) c from the 
collective expansion. 

Knowing the hydrodynamical solution as the space- 
time variations of density, velocity, and thermodynami- 
cal functions of the fluid cells, we can construct the tra- 
jectories of a pair of test pions which start initially at 
Xi and X2 and emerge finally with momenta k\ and k 2 
at freeze-out. The knowledge of the space-time trajec- 
tories of the pair of pions allows one to calculate the 
path integrals in the wave function $ of Eq. 1|12|) and 
the function P(ki,k 2 ) and P{k) after summing over all 
source elements in Eqs. © and 10- The correlation func- 
tion C(q) of the relative momentum of the two particles, 
q = \kx — k 2 |, can be constructed from P(ki,k 2 ) and 
P(ki)P(k2) by i ntegrating over the average momentum 
(ki + k 2 )/2 The HBT radius R can then be ex- 

tracted by parameterizing the correlation function C(q) 
as 



C(q) = 1 + \e- q2R2 



(17) 



In our model calculation, we shall consider the freeze- 
out configuration to be characterized by a freeze-out tem- 
perature of Tf = 0.5T . We calculate first the correlation 
function C{q) for the case (a) when the detected pions 
are assumed to originate only from the freeze-out con- 
figuration. This is the configuration considered in the 
conventional description of intensity interferometry with 



collective expansion where one assumes a chaotic source 
at freeze-out due to random pion rescattering. 

We calculate next the case (b) using the present for- 
malism to take into account the effects of multiple scat- 
tering and collective expansion with quantum probability 
amplitudes in a path-integral framework. We use the ex- 
perimental absorption cross section |l4j represented by 
the model of the absorption of a pion by a nucleon form- 
ing a delta resonance. For further comparison, we cal- 
culate C{q) for the case (c) in which the detected pions 
originate from the source under the hydrodynamical ex- 
pansion with no absorption. 

Figure 2 shows results of the two-pion interferometry 
for these cases. The conventional freeze-out configuration 
[case (a)] leads to an extracted HBT radius of 6.87 fm. 
When there is no absorption of pions [case (c)], the HBT 
radius of 4.32 fm is small. The HBT radius obtained by 
using the present formalism [case (b)] leads to an HBT 
radius of 5.31 fm which lies between the radii of case (a) 
and case (c). The extracted HBT radius depend on the 
description of the pion source and absorption. 

We conclude that for a pion source of the type con- 
sidered, the HBT radius extracted from the freeze-out 
configuration is substantially greater than the radius ob- 
tained when effects of multiple scattering and collective 
expansion on the HBT interferometry have been properly 
taken into account using quantum probability amplitudes 
in a path-integral formalism. As the latter description is 
a more realistic description of the HBT interferometry, 
the conventional assumption that the HBT interferome- 
try measures the distribution of the freeze-out configura- 
tion is therefore subject to question. 

We have carried out model calculations for some ide- 
alized configurations in order to illustrate the main fea- 
tures of the effects under consideration. These features 
will remain the same for more realistic space-time con- 
figurations. In this paper the hadronic gas is taken to be 
an ideal gas. For a more realistic case, one should con- 
sider the volume correction 0, 0] , which will lead to a 
slightly lower temperature and lower density of hadronic 
gas but will hardly affect the main features of the HBT 
interferometry considered here. The freeze-out tempera- 
ture Tf — 0.5Xb considered in this paper is lower than the 
freeze-out temperatures T f = 0.7T C and 0.9T C (T c = 160 
MeV) considered for a zero net baryon density case 
because the cross section between pion and baryon is 
greater than that between pions. The HBT radii of the 
three cases considered in this paper will increase when 
the freeze-out temperature decreases but the main fea- 
ture concerning the relative magnitudes of the HBT radii 
is unaltered. 
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